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ABSTRACT 



Context. The carbon monoxide (CO) ro- vibrational emission from discs around Herbig Ae stars and T Tauri stars with strong ultra- 
violet emissions suggests that fluorescence pumping from the ground X'Z^ to the electronic A'll state of CO should be taken into 
account in disc models. 

Aims. We wish to understand the excitation mechanism of CO ro- vibrational emission seen in Herbig Ae discs, in particular in tran- 
sitions involving highly excited rotational and vibrational levels. 

Methods. We implemented a CO model molecule that includes up to 50 rotational levels within nine vibrational levels for the ground 
and A-excited states in the radiative-photochemical code ProDiMo. We took CO collisions with hydrogen molecules (H2), hydrogen 
atoms (H), helium (He), and electrons into account. We estimated the missing collision rates using standard scaling laws and dis- 
cussed their limitations. We tested the effectiveness of UV fluorescence pumping for the population of high- vibrational levels (v=l-9, 
7=1-50) for four Herbig Ae disc models (disc mass Mdisc = 10"^, lO"'* and inner radius i?disc = l, 20 AU). We tested the effect of infrared 
(IR) pumping on the CO vibrational temperature and the rotational population in the ground vibrational level. 

Results. UV fluorescence and IR pumping impact on the population of ro-vibrational v > 1 levels. The i' = 1 rotational levels 
are populated at rotational temperatures between the radiation temperature around 4.6 /jm and the gas kinetic temperature. The UV 
pumping efficiency increases with decreasing disc mass. The consequence is that the vibrational temperatures Tyib, which measure the 
relative populations between the vibrational levels, are higher than the disc gas kinetic temperatures (suprathermal population of the 
vibrational levels). The effect is more important for low-density gases because of lower coUisional de-excitations. The UV pumping 
is more efficient for low-mass (Mdisc < 10"^ Mo) than high -mass (Mdisc > 10 ' Mq) discs. Rotational temperatures from fundamental 
transitions derived using optically thick '^CO v = 1 - lines do not reflect the gas kinetic temperature. Uncertainties in the rate 
coefficients within an order of magnitude result in variations in the CO line fluxes up to 20%. CO pure rotational levels with energies 
lower than 1000 K are populated in UTE but are sensitive to a number of vibrational levels included in the model. The '^CO pure 
rotational lines are highly optically thick for transition from levels up to £'„ppej=2000 K. The model line fluxes are comparable with 
the observed line fluxes from typical HerbigAe low- and high-mass discs. 

Key words. Stars: protoplanetary disks, molecular processes, radiation mechanisms, radiative transfer. 
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1. Introduction 

Terrestrial planets form in the inner region of protoplanetary 
discs (typically R <3-5 AU), where the gas reaches tempera- 
tures of a few hundred to a few thousand degrees Kelvin and the 
dust grains are warm enough suc h that water is not frozen onto 
their surfaces (iKamp et al.ll2010l ). The growth of solid cores in 
the inner disc is slower than in the outer disc, where grains are 
coated with an icy mantle. As a consequence, the solid bodies 
never reach the critical mass to attra ct gravitationally the disc 
gas before the gas disc has dissipated (lArmitagell2010l) . 

At a few hundred to a few thousand degrees, molecules 
are excited to their ro-vibrational levels and emit in the near- 
and mid-infrared. Probing the warm gas emissions requires high 
angular-resolution observations in the infrared. Carbon monox- 
ide is one of the most abundant species in protoplanetary discs, 
and its ro-vibrational transitions around 4.6 micron (M-band) 
are the most commonly detected lines with high spectral resolu- 
tion adaptive-optics (AO) assisted spectrometers around T Tauri 



(iNaiita et alJ 120031: iRettig et alJ 120041: ISalvk et alJ l201 ih and 
Herbig Ae star s (Blake & Boogertll2004t ICarmona et al.l 120051: 
IBrittain et al"l2009l:ivan der Plas et al.l2009l) . The CO peak emis- 
sion is also seen displaced from the source cen t re, suggesting the 
pre sence of a n inner hole ( IBrittain et al.ll2003t iGoto et al.ll2006l; 
Brown et al.l 12012) . The spectro-astrometric technique further 
improves our understanding of the emission by reaching spatial 
resolu tions of the order of 0.1-0.5 AU for nearby (100-140 pc) 
discs (iPontoppidan et al.ll2008l) . iGoto et all (|2012|) observed the 
CO ro-vibrational emission from the HD 100546 disc with a spa- 
tial resolution of 10 AU at 100 pc. They found warm CO gas 
(400-500) out to a distance of 50 AU, which is evidence for the 
presence of a warm molecular disc atmosphere. Very dense and 
hot gas located within a few tenths of AU from the stars can 
be detected by the CO overtone em issions jTatuUi et al.ll2008l: 
iBerthoud et aLll2007l: iThi et al.ll2005l) . 

Line profiles are diverse, from broad single peaked to clear 
double peaked: the latter is a si gnature of Kepler i an rotation 
of the gas between 1 and 5 AU dSalvket afll2009l) . iBast et all 
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(l201lh focused on 8 out of 50 discs that show broad single- 
peaked CO profiles. They proposed that those 8 discs have either 
highly turbulent inner disc gas and/or that a slowly moving disc 
wind contributes to the emission. 

The M-band observations show lines emitted from a large 
range of excitation levels (from fundamental v = 1 ^ P(l) 
and R{Q) lines to "hot" lines (v' > 1, Av = 1) which make 
it possible to probe disc regions with different excitation con- 
ditions. CO vibrational diagrams of Herbig Ae discs indicate 
vibrational gas temperatures higher than the rotational ones (a 
few thousand K instead of a few hundred K). This suggests that 
the population of high vibrational levels (v >1) may be domi- 
nated by IR and UV fluorescence excitation dBrittain et al.l2007t 
iBrown et alJl2012h . The V > 1 ro-vibrational levels can also be 
efficiently populated by UV fluorescence for T Tauri discs, es- 
pecia lly for strongly accreting systems with excess UV emis- 
sion (iBast et alJl201 ih . Fluorescence excitation in T Tauri discs 
is supported by Hubb le Space Telescope observations of UV- 
pumped CO emission dFrance et alj|201 it) . 

Many models of CO ro-vibrational emission made the as- 
sumpti on of local thermodynamical equilibrium (LTE) popu- 
lation dHugelmever et al J 120091: iRegalv et al.l 120101: ISalvk etalJ 
l20TTI) . Nevertheless, a detailed understanding of the CO ro- 
vibrational lines, especially the spatial location of the lines and 
the efficiency of the IR/UV fluorescence, requires non local 
thermodynamic equilibrium (NLTE) modelling and a chemo- 
physical code that includes detailed continuum and line radia- 
tive transfer. The accuracy of NLTE modelling depends on many 
factors, including the availability of accurate collision rates. The 
difficulties associated with the measurement and computation of 
accurate rates result in sparse trustworthy data. Completeness 
in the rate coefficients can be attained only by using approxi- 
mate extrapol a tion ru les derived from physical considerations. 
iKrotkov et ai] ( Il980l) and IScoville et all ( Il980l) studied the ef- 
ficiency of UV and IR pumping but considered only pure- 
vibrational levels. 

In this paper, we describe the implementation of a large CO 
model molecule, which includes many ro-vibrational levels in 
the ground electronic state X and the electronic state A, into the 
code ProDiMo. The radiative chemo-physical code ProDiMo 
solves the continuum radiative transfer, the gas heating and cool- 
ing balance, the gas chemi stry, and the disc v ertical hydrostatic 
structure self-consistently dWoitke et al.ll2009l paper I). This pa- 
per is the fourth in a series dealing with aspects in pro toplane- 
tary disc modelUng dKamp et al.ll20IOl : iThi et al.ll20I iL paper II 
& III). 

We will use the code to model two typical discs around an 
A-type star: a solar-nebula-mass disc (Mdisc - 0.01 M©), and a 
low-mass disc (Mjisc - 10""* M©). The choice of these models 
isq justified by the expectation that the strength of UV/IR flu- 
orescence will depend on the density and optical depth in the 
discs. 

The paper is organized as follows: we present and discuss 
the CO spectroscopic and collisional data implemented in the 
ProDiMo code in Sec.|2l The main feature of the ProDiMo code 
and the disc model parameters are described in Sec. l3.1l and l3.2l 
In Sec. 13.31 we discuss the model outputs. Finally, we present 
our conclusions and recommendations on the interpretations of 
CO ro-vibrational observations in discs in Sec.|4] 



2. CO molecular data 

2. 1 . Spectroscopic data 

iTashkun et al.l d2010l) compiled and evaluated existing experi- 
mental frequencies for transitions within the ground electronic 
state. Their work i s based on accurate experimental frequenc y 
measurements (e.g. lGeorge et alj|l994l : lGendriesch et al.ll2009h . 

We compared their ene rgy levels and frequencies with those 
from lChandra et all d 1 9961) and lHure & Roueffldl993h an d found 
no significant differences. In addition. IChandra et al.l dl996h 
computed Einstein spontaneous probabilities A'll for all ro- 
vibrational transitions between 140 rotational levels and 10 vi- 
brational levels in the ground electronic state. We adopted their 
trans ition probabilities, which compared well with other studies 
(e.g. lOkadae"tani2002h . The transitions between the electronic 
excited ro-vibrational levels A and the ground ro-vibrational lev- 
els Z^S^ (4th positive system), which occur around ~ 1600 A, 
were derived from the b and-averaged transition probabilities of 
Beegle et~al1 dl999h and lBorges et al.1 d2001h using t he formula 
and Holn-London factors in iMorton & Noreaul dl994). which are 
given in Table [T] Figure [T| shows a sketch of the fluorescence 
mechanism involving a ro-vibrational level in the ground elec- 
tronic state X'S^ and a ro-vibrational level in the excited state 

A'n. 



Table 1. Holn-London factors taken from iMorton & Noreaul 
(il994t) . 
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2.2. Collisional data 

Although CO is one of the most studied interstellar molecules 
with a large amount of experimental and theoretical collision 
rates available, the lack of completeness of any individual study 
obliged us to use scaling laws to fill the gaps. Since protoplan- 
etary discs show a large range of physical and chemical proper- 
ties, we need to consider collisions of CO with H, He, H2, and 
electrons. In the rest of the section, we will first review existing 
theoretical and experimental collision rates before discussing the 
strengths and limitations of scaling laws. 

2.2.1. CO-H collision rates 

The v=l— >0 CO-H de-excitation rate coefficients were first 
measured by iMillikan & White! (Il963h in shock-tube exper- 
iments. Subsequent works include the shock-tube experi- 
ments by Ivon Rosenberg etlD (ll97lh.lGlass & Kirondg dl982h . 
and more recentlv IkozIov et aiT l 2000h (between 2000 and 
3000 K). Shock-tube experimental measurements are de- 
excitation timescales,which can be fitted in terms of the Landau- 



Teller rate coefficients in cm s ' dAvres & Wiedemannlll989h : 



= 4.2 X 10 



_i9exp[B-0.0069A/3i/3] 



(1) 



/3{l-e-l^) 

with P = Ey^i/T and the val ues for A and B given in Table |2] 
iNeufeld & Hollenbachldl994 gave a simple formula that fits the 
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Fig. 1. CO energy potential curves on the left and a schematic of the UV fluorescence mechanism on the right. We are concerned 
with transitions between the ground electronic state X'S^ and the electronic excited state A '11. The arrows show a possible UV- 
fluorescence pumping from a ground electronic level (X'S^) at v" = 1 to a level in the electron ic state A'll, followed b y a de- 
excitation to the level v" = 2 in the ground electronic state. The value for the curves are taken from lCooper & Kirbvl (Il987h . 



shock-tube data for T > 300 K: 

ky{\ -> 0) = 7.9 X 10^'^rexp (-12087). 



(2) 



This formula supersedes an older one given by 
iHoUenbach O^cK ee ( 1979): 



/t,(1^0) = 3xlO-'2Vfe-<2x'°W 



-3080/r 



(3) 

Experimental values for transitions oth er than the v=l— >0 are 
not available. iBalakrishnan et al] ( l2002h performed quantal scat- 



Table 2. Landau-Teller coefficients for CO-partner collisions de- 
rived from shock- tube experiments. 
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87 
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Notes, d) iGlass & Kirondl (Il982h.(2) iMillikan & Whiia jl963h. G ) 
Ivon Rosenberg et alj ( ll971h : lKozlov et all ( l200(3h . (4) lMillikanl ( 119641) . 



tering computations for rotational and band-averaged vibrational 
collision rates between CO and hydrogen atoms. They used the 
close-coupling method for the rotational transitions and the in- 
finite order sudden approximation (lOS. iFioweJ l2007l) for the 
vibrational transitions for a large range of gas temperatures (5 
< T < 3000 K). The lOS approximation reduces the computa- 
tional effort but is strictly valid only for energetic collisions, i.e. 
AE <& kT, where AE is the energy difference between two lev- 
els. The close-coupling computations were done for rotational 
levels up to 7=7. The vibrational transition rates were computed 
up to V = 4. 

Figure|2]shows the different rate coefficients. The most strik- 
ing feature is the lar ge dispersion (a few orders of magnitude) 
between the sources. 



Avres & Wiedemamil (Il989l) chose the val 

ues of iGlass & K ironde' ('1982'), which a re mu ch larger than 
the older values of von Rosen berg et al.Jl97lh . However. 

.1 (I2000h 



cent measurements byLK ozlov et al 



re- 



reproduced the val- 



ues of Ivon Rosenberg et alj (|197T1) . Interestinglv. the theoretical 



10-' 
10-^" 

i 10-" 



o Rosenberg - Kozlov (2000-3000K) 
m Miliikan & White 
+ Glass & Kironde 




Baial<rishnan 
Hoilenbach & IVlcKee 
Neufeld & Hoilenbach 



500 



1000 



1500 
T(K) 



2000 2500 



3000 



■*0) rate coefficients of 
The large discrepancies 



Fig. 2. Vibrational de-excitation (v=l 
CO by collisions with H (in cm-' s '). 
between the experimental values may stem from the reactiv- 
ity of H-atoms with CO. We adopted the theoretical rates from 
IBalakrishnan et all ( 2002h. which match the experimental values 
of IGlass & Kirondei (119821) (plus signs). 



values calculate d by Balakrishnari et al. (l2002l) are much closer 
to the values of iGlass & Kironde' (1982) and are best ap prox- 
imated by the formula of .Neufeld & HollenbachI (|1994|) . One 
possible reason for the large discrepancies is the reactivity of 
CO with atomic hydrogen (an open shell species), which can 
fal sify the measurements. W e chose to use the theoretical rates 
bv IBalakrishnan et all (l2002h . 



2.2.2. CO-He collision rates 



ICecchi-Pestellini et al.l (l2002l) computed rotational and vibra- 
tional rate de-excitation coefficients of CO in collision with He. 
The rotation rate coefficients in the ground vibrational level are 
computed by the quantal close-coupling method. The vibrational 
rates are co mputed for 500 < T < 1300 K using the lOS ap- 
proximation. iKremsl (l2002l) calculated quantum-mechanically 
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Fig. 3. Vibrational de-excitation (v= 1— >0) rate coefficients of CO 
by colli sions with He (in cm^ s '). We adopted in this study the 
rates of iKrems] ( 2002h at temperatures < 1500 K (plus signs) and 
those of Cecchi-Pestellini et alJ (l2002l) at very high temperatures 
(> 1500 K, diamonds). 
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Fig. 4. Vibrational de-excitation ( v= 1 — >0) rate coefficients of CO 
by collisions with H2 (in cm-' s"'). Normal H2 is a mix of 0-H2 
and P-H2 with an ortho-to-para ratio at LTE. 



ro-vibrational rate coefficients fo r 35 < T < 1500 K. Rates 
were computed by [Flowed (1201 2| ) and Reid efal] d 19951) . The 
experi mental fit coefficients of Millikani (119641) are given in 



Tablel2] lwickham- Jones et aP (11987 ) reported experimental data 
from the PhD thesis of Chenery (1984), who measured rates 
down to 35 K. The experimental rate coefficients are compared 
to the theoretical values in Fig. [3] The theoretical values of 
ICecchi-Pestellini et al. I (l200l are a factor ~ 2 smaller than the 
experi mental values, whereas the values computed by iKremsl 
(12002 ') match the shock-tube data down to r= 300 K. The for- 
mula [1] that fits the experimental shock-tube data under-predicts 
the values for T < 300 K. We adopte d the values of Krems 



(|2002|) up to r= 1500 K and those of ICecchi-Pestelhni et al 
(2002) above r= 1500 K. The rate coefficients for CO in col- 



lisions with He are around three orders of magnitude smaller 
than the collision rates with H. This behaviour can be explained 
by the fact that He is a close-shell species. 

2.2.3. CO-H2 collision rates 

In the case of collisional excitation of CO by H2, we need to 
consider the colli sions with ortho-H2 and para-H2 separately. 
lYang et al] (1201 Oh computed with the close-coupling method 
rates for temperatures between 1 and 3000 K involving ground 
vibrational CO transitions and with rotational levels up to 7=40. 
Rates for vibrational transitions are scarcer and are given aver- 



aged o ver the rotational levels. Between 30 and 70 K. lReid et alj 
(Il997h showed experimental and theoretical rates for the CO 
v=1^0 transition. Between 70 and 290 K, rates derived from 
shock-tube experiments are available ( Hooker & Millikanll 19631: 
[Andrews & Simpsonl I1975L 119761: iGlass & Kirondd Il982h . At 
temperatures higher than 300 K, we can assume that the Landau- 
Teller law i s vahd and expand the sh o ck-tube rates to a few 
thousa nd K dAvres & Wiedemannll989l) . lNeufeld & HoUenbachl 
d 19941) provided a convenient analytical expression of the shock- 
tube data (r >290 K): 



k,{\ ^ 0) = 4.3 X 10-"*rexp(-68r^'^) 



(4) 



Alternatively, one can use the Landau-Teller formula with the 
values in Table |2l 



In Fig. m we compare the two shock-tube experimental 
rate coefficients, the experimental and theoretical values of 
iReid et all (Il997h . and our adopted fit 



k,{\ 



0) = 5 X 10"'''rexp(2 X lo^^r). 



(5) 



for r < 300 K and the fit to the 



bv lNeufeld & HoUenbachl (1 1994 above 300 K. We notice a large 



Hooker & Millikani (IT963I) data 



discrepancy between the experimental and theoretical rates at 
low temperatures. The shock-tube data are consistent with each 
other. It is known that the rate extrapolations to lower tempera- 
tures with the Landau-Teller law may not be correct. Collisions 
with H are around two orders of magnitude more effective than 
collisions with H2 at all temperatures. 



2.2.4. CO-electron collision rates 

The collision excitation and de-excitation rates of CO by elec- 
tr ons for rotationa l trans itions are computed using the theory 
of [Dickinson et al.l (Il977|) based on the Bo rn approximation for 
A7 = ± 1 (see also iRandell et all ri996h . We adopted a ro- 
tational constant Bco = 1.922529 cm"', a centrifugal distor- 
tion Deo - 6.1210779 X 10~^ cm ', and a dipole moment of 
0. 1 12 Debye. iRistic et all d2007l) measured the vibrational rates 
of CO by collisions with electro ns for vibrational levels up to 
ten. iMorgan & TennvsonI (Il993l) also comp uted cross sections 
for vibrational transitions. iThompsoiil dl973h gave a simple for- 
mula for the v= 1 — »0 rate 



yt,,(l ^ 0) = 1.4 X 10"V'^^(l 19e"^-^2^(l -H 4.22/3)), 

(6) 

where again ji — Ey^\IT . An alternative formula is given by 
iDraine&Robergdd 19841) 



k,(v =1^0) = 1.9 X 10"" Vf 
10 



2420 



, 10500 / 10210\) 
1 H exp 



(7) 
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Fig. 5. Vibrational de-excitation {v- 1 — >0) rate coefficients of CO 
by collisions with electrons (in cm^ s '). 




10000 
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Fig. 6. Extrapolated CO - 0-H2 pure rotational rate coefficients. 
The X-axis corresponds to the upper level energy in Kelvin. The 
diamonds are rates from lYang et al.l (l2010h . 



Finally, we propose our ow n analytical formula, which fits the 
values of lRistic et aP (l2007h 



A:,(v= 1 ->0)= 1.9xl0""r'/^ 
7 2000 \'^^^ 



law for the base rates) formula: 



50. 1 



70000 / 10210\ 
exp ' ' 



— ; 



exp(-4.5 X lO^^r). (8) 



The measurements for the v=1^0 transition are shown in Fig.|5] 
and compared to the three analytical fitting formulae. The 
iThompsonI d 19731) formula deviates significantly from the ex- 
perimental data. We used a combination of formula (O, which 
fits the experimental data relatively well down to 400 K.Below 
400 K, we adopted a flat rate coefficient. 



2.3. Rafe coefficient scaiing iaws 

The current experimental and theoretical works do not provide 
all the rate coefficients needed to model the CO statistical ro- 
vibrational population of the ground and excited A ' n elec- 
tronic levels. Therefore, we have to resort to extrapolation and 
scaling laws to fill the gaps. 



2.3.1. Extrapolating rotational transition rates 

A few extrapolating laws attempt to estimate rates from ex- 
isting ones. We used the energy-corrected sudden (ECS) scal- 
ing law wit h a hybrid exponential power-law for the bas e 
(kj-.^o) rates (iDepristo et al.lll979l: lDepristol[T97l lGreenlll993h . 



iGoldflam et al.l d 1977 ) showed that rates between any levels can 
be derived knowing the rates to the ground level within the lOS 
approximation. To overcome the limitatio n of colhsion energies 
greater than the level energy differences, IDepristo etal.l(ll979h 
proposed a correction for the inelasticity and finite collision du- 
ration effects using second-order perturbation theory for low- 
collision energies (i.e. for collisions at low temperatures). A 
de-excitation (/' > J") rate is given by the ESC-EP (energy- 
corrected sudden scaling law with a hybrid exponential power- 



(2 J" -H l)exp 



(Ej. - E,r ) 



kT 



YP-L + mi} ^ Q Q q' j ^o,y'^o,o, (9) 



where 



J' L J" 




(10) 



is the Wigner 3- j symbol and Ejt and Zsy^y-i are the upper and 
lo wer level energy, respe ctively. The correction factor introduced 
bv IDepristo et al.l d 19791) is defined by 



I + Trie' 



with 



1/2 



(11) 



(12) 



where yu is the reduced mass in atomic mass units, is the 
inelasticity to the next lowest level {Ejt - Ej»^j>-i) in cm"', 
and Ic is an adjustable scaling length in A. We adopted 1^=3 A 
(ISchoier et al.ll2005l) . The ESC formula is recovered for Ic-O A. 

The scaling law requires the knowledge of the base rates 
[J — > 0). Collision rates with hydrogen atoms to the ground level 
are known for J up to 7, with H2 up to 40, and only for 7=1— >0 
for electrons. Further estimates of the rates to the ground level 
are needed. The ESC-E P model provides an appro ximation for 
rates to the ground level (iGreenll 1993t lMilio3[T990h : 

^o,y'^o,o = A{T) [J'iJ' + l)]-^exp(-/3Ej,/kT), (13) 

where A(T), y, and fi are free parameters, which are estimated by 
fitting equation[T3]to the existing data using a least-square solver 
for over-constrained systems. We performed the extrapolation 
for collisions with 0-H2, p-Hi, He, and H but not with electrons. 
An example of extrapolated rate coefficients is shown in Fig. |6] 
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2.3.2. Extrapolating vibrational transition rates 

The extrapolation of an existing rate to low temperature was per- 
formed by fitting the available rates to the formula 



ky - ax exp(r + b). 



(14) 



where a and b are the two fitting parameters. 

The extrapolation of rates between two vibrational levels 
without any known rate follows physically motivated scaling 
laws. Based on the initial work of iProcaccia & Levind (1 19751) . 
IChandra & Sharmal (1200 ih proposed using the Landau-Teller re- 
lationship for transitions between adjacent states (Av =1) 



k,{v' v' - 1) = v'Kil 0), 



and for Av > 1 : 

k,{v' 



V") = T^^'Xil 



2(v' -2) + l 
lEUtzud (Il983h argued for a different formulation 



0) 



(15) 



(16) 



k,{v' v") = {v'-\'")k,il 0) 

exp(-(v'-v"-l) 



\.5E,lkT 
1 + 1.5E,/kT 



For Av - 1 (v" = v' - 1), this last formula reduces to 
k„iv' v' - 1) = k,{l ^ 0). 



(17) 



(18) 



IScoville et alj ( Il980l) argued that the probabilities of 
colUsion-induced vibrational transitions are proportional to 
the corresponding radiative transition matrix elements (Born- 
Coulomb long range interaction). Therefore, the vibrational de- 
excitation rate coefficients from v' to v" can be derived from the 
V = 1 — > rate coefficients: 



k,(v' 



v") = ^kAl 



0), 



(19) 



Mo 



where A,,/,," is the vibrational band transition probabilities. The 
three formulations give quite different rate estimates. In order to 
choose which one is the most appropriate, we compared the pre- 
dictions with detailed calculations and/or experimental values. 

Figure Q shows the rate coefficient ratios for transitions be- 
tween adjacent levels in c ollisions with atom ic hydrogen using 
the detailed calculation of lBalakrishnan et"aD (|20Q2j. The ratios 
are grouped around one with a small spread. The rate coeffi- 
cient r atios suggest that the expansion rules proposed bv lElitzurl 
(Il983h are more appropriate to expand the detailed rate coeffi- 
cients. 

The rate coefficient ratios for CO-He collisions vary dramat- 
ically with temperature (Fig.[8]i. For temperatures below 1000 K, 
the ratios are close to 2 x v ' , a va lue that corresponds nei- 
ther to the IChandra & Sharma ( 200 lb prescription (= v') nor to 
that championed bv ^Elitzuj ( 1983h (=1). We chose to use the 
[Chandr a & Sharma (2001) for Av >1 and the factor 2 x v' for 
Av=l. 

The choice is open concerning the rules to apply to estimate 
the rate coefficients between CO and H2 and between CO and 
electrons other than \'"=1— >0. We chose to use the same combi- 
nation of scaling rules as for CO-He collisions. We estimate the 
scaling rules to provide rates correct within an order of magni- 
tude. 




200 400 



600 
T(K) 



800 1000 1200 



Fig. 7. Vibrational de-excitation rate coefficient ratios k(v'^w'- 
1)/ ^(1— >0) of CO by H . Rate coefficients are taken from 
iBalakrishnan et al.l ( l2002l) . 
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Fig. 8. Vibrational de-excitation rate coefficient ratios k( v'—^v'- 
1)/^(1 — >0) of CO by He. Rate coefficients are taken from lKremsl 
(l200l . 



2.3.3. Estimating ro-vibrational transition rates 

The vibrational rates are averaged over all ro-vibrational tran- 
sitions. To obtain individual rates between ro-vibrational states, 
we apply the method described in lFaure & Josse lin (2008'). This 
method assumes a complete decoupling of rotational and vibra- 
tional levels. The state-to-state ro-vibrational rates are related to 
the corresponding rates for the pure rotational transitions in the 
ground vibrational state: 



J'^v",J" 



Pv'v"ko,J'^0,J" 



The factor Pv'v" is defined as 



P.' AT) = 



■IjJ' 8J'^Wi-Ev',j'/kT) 



(20) 



(21) 



This formula ensures detailed balance. 



2.3.4. Approximating electronic transition collision rates 

In the absence of collision rates connecting the ground X'E^ to 
the A'n electronic levels, we adopted the rough values for CO- 
electron collision rates provided by the g approximation. The g 
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(Ivan Regemorteil[l962h approximation relies on the Born-dipole 
(also called Born-Coulomb) long-range interaction approxima- 
tion and relates the collision strength to the transition pro babil- 
ity and wavelength for allowed transitions (Iltikawall2007h . The 
Born-dipole approximation is valid for species with a permanent 
dipole, which is the case for CO {fi-Q.l22 Debye). The collision 
rate for allowed transitions is given by 

8.6291 x 10-''^ 
kui = p r (22) 

For neutrals, the g-bar value can be approximated by 

-g - (kT/hv) (23) 

The collision strength F is related to g-bar by 

F = 2.388 X 10-'' Al^gA.ig (24) 

The approximation provides rate coefficients correct within an 
order of magnitude. 

2.3.5. Combining tlieoretical, experimental, and estimated 
rates 

The resulting CO rate matrix after combining all existing and es- 
timated rates is complete in the sense that all ro-vibrational lev- 
els in the ground electronic state connected by a radiative transi- 
tion are also connected by collisions. Therefore all levels in the 
ground electronic levels can be populated in LTE at very high 
densities and temperatures. We have neglected the radiative and 
collisional transitions between levels in excited electronic states 
since radiative transitions to the ground electronic state are very 
fast. The "completeness" is required to ensure that we don't arti- 
ficially create sub-thermally populated levels. The methodology 
described above was implemented in the radiative photochemi- 
cal code ProDiMo. The main collisional partner in ro-vibrational 
transitions is the atomic hydrogen by a few orders of magnitude. 



Table 3. Disc and C O parameters. The parameters are defined in 
IWoitkeetal.l(l2009h 



stellar mass 


M, 


2.2 Mo 


stellar luminosity 


L. 


32 Lo 


effective temperature 




8600 K 


distance 


d 


140 pc 


disc inclination 


i 


° (face-on) 


total disc mass 




10-2, 10-* Mo 


disc inner radius 


^in 


1,20 AU 


disc outer radius 


^out 


300 AU 


vert, column density index 


e 


1 


inner rim soft edge 




on 


gas to dust mass ratio 


5 


100 


dust grain material mass density 


Pdusl 


3.5 g cm-^ 


minimum dust particle size 


^min 


U.Uj /jm 


maximum dust particle size 




1000 yum 


composition 




ISM silicate 


dust size distribution power-law 


P 


3.5 


Ht cosmic ray ionization rate 


^TcR 


1.7 X 10-" S-' 


ISM UV field w.r.t. Drame field 


X 


1.0 


abundance of PAHs relative to ISM 


/pah 


0.1 


a viscosity parameter 


a 


0.0 


lUlUUlCIlCC WiUlll 


6v 


0.15 km S-' 


reference scale height 


Ho 


15 AU 


reference radius 


Ro 


100 AU 


flaring index 


P 


1.2 


COX'S+ 


rot. levels 50 


COA'n 


rot. levels 50 


COX'E+ 


vib. levels up to 9 


COA'n 


vib. levels up to 9 



all the discs. The inner radius was set to two values to model 
the presence of an inner disc hole and its effects on the CO ro- 
vibrational emissions. 



3. CO fundamental and hot-band emission from 
l-ierbig Ae discs 

3.1. The ProDiMo code 

A det ailed description of the code can be found in (IWoitke et al.l 
20091). Subsequen t additions to the cod e are documented i n 
Kamp et al.l fcOlOl) . iThi et aH (l201 ih . and IWoitke et~an (1201 ih . 
Although not used in this study X -ray heating and X -ray chem- 
istry are included in ProDiMo (lAresu et al.l 1201 ih . We deal 
in this paper with the excitation mechanisms of the CO ro- 
vibrational levels and do not attempt to fit actual observations. 

3.2. Disc model description 

We chose to model two representative discs, which can be 
found around Herbig Ae stars. Our generic star has an effec- 
tive temperature of 8600 K, a mass of 2.2 Mq, and luminosity 
of 32 Lq. We have no intention to fit a particular object in this 
study. The input stellar sp ectrum is t aken from the PHOENIX 
database of stellar spectra dBrott & H auschildt 2005). The disc 
model parameters are summarized in Table[3] The disc geometry 
is described by a single zone power-law surface-density profile 
with index e, an inner radius Ri^, an outer radius Rout, a disc gas 
mass Mdisc, and a standard gas-to-dust mass ratio A of 100 for 



The dust grains are de fined in both discs by their composition 
(astronomical silicates. iLaor & Draind l 1993h and their size dis- 
tribution (a power-law ranging from flmin to Umax)- We did not it- 
erate on the vertical hydrostatic profile but used instead a power- 
law H = //()(/?//?o)-'^ to model the gas scale-height H. We chose 
to model two flaring discs (flaring index p=l.2) because obser- 
vations suggest that the CO levels populated by UV-fluorescence 
pumping seem to occur in flaring discs. 

The two discs differ from each other only in the total 
(gas+dust) disc mass Mdisc- The first disc has a mass consistent 
with the typical solar nebula mass of 0.01 Mq, while the sec- 
ond disc is more akin to an object in transition from a primordial 
gas-rich disc to a secondary debris disc with a total mass of 10-^ 
Mq. The choice of the two masses is motivated by the possibil- 
ity to assess the effects of infrared and UV fluorescence on the 
CO population in discs with vastly different density structures. 
The discs are assumed passively heated only (i.e. without vis- 
cous heating: a - 0), and we do not consider X-ray heating be- 
cause the X-ray luminosity is generally low compared to the UV 
luminosity in Herbig Ae stars. X-ray heating and c hemistry can 
play an important role in discs aro und T Tauri stars ( lAresu et al.l 
1201 lll2012l:[Meiierink et al.ll2012l) . 

The main source of heating for the gas is the stellar UV 
via photoelectric effects on polycyclic aromatic hydrocarbons 
(PAHs) assumed here to be circumcoronenes (C54H18). The 
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abundance of PAHs is assumed at 10% of the standard inter- 
stellar abundance (/pah - 0.1). 

The gas and dust were assumed well-mixed with no dust 
settling. A detailed discussion on the heating an d coohng in 
HerbigAe discs can be found in lKamp et al.l (|2010|) . 

The continuum radiative transfer and dust thermal balance 
are calculated in two dimensions usi ng th e long-characteristic 
method ( Woitke et al. 2009; Thi et a l...201 lb . The radiative trans- 
fer numerical impleme ntation was benc hmarked against other 
radiative transfer codes dPinte et al.ll2009h . The gas temperatures 
were computed by balancing the heating and cooling (mostly 
radiative) processes. The CO self-shielding against photodisso- 
ciation uses factor s interpolated from values in look-up tables 
dVisser et al.ll2009l) . 




10^ 10' 10° 10' 10' 10^ 10" 10' 
optical depth 

Fig. 9. Escape probability for a semi-infinite slab. We adopted 
a - 10^^ for illustration. 



CO ro-vibrational transitions have moderate Einstein prob- 
abilities (A =1-10 s ') so that the line transfer occurs in the 
Doppler core of the line profile. On the other hand, UV- 
fluorescence pumping occurs via electronic transitions (A'll- 
X'E+) with high Einstein probabiHties (A^IO'^-IO^ s"') and thus 
with large natural line width. The line core becomes quickly op- 
tically thick (line self-shielding) with the UV line transfer hap- 
pening in the Lorentzian wings of the line. Therefore the use of 
the Voigt profile in the line transfer is warranted. 

The CO levels are computed using a 1 + lD escape proba- 
biUty method for Voigt line profil es. We adapted the analyti- 
cal formulation of Apruzesd (1 19851) to match the formulation in 
IWoitke et al.l (l2009l) at low optical depths and small intrinsic line 
widths. At high optical depth, the escape probability /3 varies as 
T V2 (Voigt profile) instead of r ' (Gaussian profile): 



[0.5(1 + 1.5r)-i if T< 1 
/3 = I 0.25r-i i5 if 1 < T < 

0.25T-° Vr" ''5 if r > Tr, 



(25) 



where a is defined as r/(47rAv£)) and r^. is the critical optical 
depth and is defined as 0.83/(fl(l -i- V^))- The sum of the natural 
and collisional width is F while the effective Doppler width is 



Avd (iRvbicki & Lightmanl[T986l) . The escape probability func- 
tion for a Voigt line profile is illustrated in Fig. |9]for the case 
a = 10"^^. In addition to the UV line shielding, dust grains con- 
tribute strongly to the UV flux attenuation. The main limitation 
in using the escape probability technique is that it does not take 
overlapping line effects into account. 

Excited vibrational levels can also be populated by the ab- 
sorption of an IR photon emitted by the star and the dust grains, 
a phenomenon called IR pumping. We also checked the effect of 
IR pumping on the vibrational ground rotational population by 
modelling two discs with /?in=l AU, one with Mdisc = 10"~ Mq 
(models Ic, Id, and le) and with Mdisc=10""* Mq (models 3c, 
3d, and 3e). For each disc model, we computed the rotational 
level population by taking into account: (1) ro-vibrational levels 
up to v=5 (models Ic and 3c), (2) ro-vibrational levels with v=0 
and v=l (models Id and 3d), and (3) the rotational levels in the 
ground vibrational only (models le and 3e). 

In total, we ran four series of models with and without the 
electronic levels and with different numbers of vibrational levels 
to assess the importance of UV fluorescence and IR pumping 
(see Table |4]i. Complete spectra from four to five microns are 
generated for discs seen face-on (zero inclination). 

3.3. Results & Discussion 

3.3.1. CO chemistry and location of tine CO ro-vibrational 
emission 

We focus on the CO chemistry of model la. The disc density, 
dust temperature, and gas temperature structures are shown in 
Fig. [To] The panels in Fig. [TTIshow the enhancement (x) with 
respect to the interstellar UV, the abundance relative to the total 
number of hydrogen-nuclei for He, electron, atomic hydrogen, 
molecular hydrogen, vibrationally excited H2 (H2exc), C^, C, 
CO, CH4, OH, and H2O. The transition from C+ to C occurs at 
logOif/n) between -1 and -2, and the transition between C and CO 
occurs at \og{x/n)~ -2.5 to -3. 

In Fig. [121 we show the location of the CO \' = 1 - P(19) 
and V - 4-3 P{20) emissions overlaid with the number density 
of the three main collision partners (H, H2, and electrons) for 
a 10"~ Mq disc model taking UV pumping or not into account 
(models la and lb). Similar plots are displayed in Fig. [13] for 
models 3 a and 3b. The CO lines are emitted up to a few AU after 
Rin, where the density is high enough such that the level popu- 
lation is the closest to the LTE population for the fundamental 
transitions. The emitting area becomes more and more confined 
to the inner rim as the upper energy of the transition increases. 
In the 10"" M0 disc model the fluxes are a few percents higher 
when UV pumping is switched on. The fundamental P(19) line 
is emitted in region where the gas is between 500 and lOOOK. 
In the vertical direction, the CO ro-vibrational lines are mostly 
emitted in the disc region with z/r~Q.2 and R< 5 AU. 

In the emitting region, the extinction in the visible Ay is be- 
low 1 and some H2 molecules are excited by UV and/or hot 
enough to be vibrationally excited (the excited H2 abundance 
reaches lO^'^-lO""^). As soon as it is formed, the warm or vibra- 
tionally excited H2 reacts quickly with atomic oxygen to form 
OH. In turn, OH reacts with H2 to form H2O or with atomic 
carbon to form CO. 

The OH and water vapour abundance are high for z/r~Q.2 
and R< 5 AU. The central role playe d by O H in hot gas chem- 
istry is discussed in model details in lThi & Bik (2005). The la 
and lb disc models are massive enough such that in the inner 
disc midplane [z/r <0. 1) the most abundant carbon-bearing gas- 
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Table 4. CO LTE rotational temperature and vibrational temperature for the UV/no-UV models derived from analyzing the CO 
fundamental emission rotational diagrams, max(rdust) is the maximum dust temperature attained in the disc. n.a. means not available. 



Model Md,sc Rn, A'n UV T,™ T,,t, max(rd„„) 

number of number of pumping? v = 1 — > v=0,...,5 





\1\IQ } 


I'M 11 




vih IpvpIq 










la 


10-2 


1 


9 


9 


yes 


111 


1413 


857 


lb 


10-2 


1 





9 


no 


111 


941-1413 


857 


Ic 


10 


1 





5 


no 


n.a. 


1413 


857 


Id 


10-2 


1 





2 


no 


n.a. 


n.a. 


857 


le 


10-2 


1 





1 


no 


n.a. 


n.a. 


857 


2a 


10-2 


20 


9 


9 


VPS 


655 


2025 


184 


2b 


10-2 


20 





9 


no 


655 


576 


184 


3a 


10-" 


1 


9 


9 


yes 


907 


1740 


815 


3b 


10-" 


1 





9 


no 


907 


858-1192 


815 


3c 


10-" 


1 





5 


no 


n.a. 


858 


815 


3d 


10-" 


1 





2 


no 


n.a. 


n.a. 


815 


3e 


10-" 


1 





1 


no 


n.a. 


n.a. 


815 


4a 


10-" 


20 


9 


9 


yes 


658 


3177 


172 


4b 


10-" 


20 





9 


no 


658 


560 


172 


5a 


10-3 


1 


9 


9 


yes 


750 


1438 


838 


5b 


10-3 


1 





9 


no 


750 


1096 


838 


6a 


10-3 


5 


9 


9 


yes 


772 


1731 


352 


6b 


10-3 


5 





9 


no 


772 


1000 


352 


7a 


2.5 X 10-3 


5 


9 


9 


yes 


734 


1628 


357 


7b 


2.5 X 10-3 


5 





9 


no 


734 


915 


357 


8a 


5 X 10-3 


5 


9 


9 


yes 


726 


1575 


360 


8b 


5 X 10-3 


5 





9 


no 


726 


872 


360 




r[AU| r|AU| r[AU] 

Fig. 10. Density, Tdust, and Tgas structures for model la. The contours in black dashed lines in the density panel show the location 
of Av=l and 10 from z - oo. The red dashed lines mark the location of the radial Ay^l calculated from T^in. In the middle panel 
various dust temperature contours are shown in black dashed lines. In our models, Tdust and Tgas are not equal. 



phase species are methane and C2H2 {R< 3 AU), whereas most 
of the oxygen is locked in water vapour up to 5 AU. The rem- 
nant oxygen is the form of atomic oxygen, CO, and other minor 
species. The gas in the inner disc is dense enough such that UV 
fluorescence is quenched by collisional de-excitations between 
the ro- vibrational levels. 

3.3.2. CO Spectral Line Energy Distribution 

The analysis of the CO ro-vibrational lines can be performed ei- 
ther by comparing directly the observed and modelled line fluxes 
or by using rotational diagrams. In Fig. fT4l and [T6l we plotted 
the continuum-subtracted line fluxes for model 1 to 4 (spectral 
line energy distribution, SLED) for the fundamental (v=l-0) and 



one "hot" transition (v=4-3). The shape of the SLEDs varies 
from model to model and between models with and without UV 
pumping. The large variety of the SLEDs illustrates the sensi- 
tivity of the CO ro-vibrational line fluxes to the disc parameters 
(A^disc, ^in, ctc ...). The line fluxes range from IQ-''' W mr~ to 
a few IQ-'^ W mT^ for a low-mass disc (IQ-" Mq) with an in- 
ner hole to a few IQ-'* W mT^ to for a massive disc (IQ-^ M©) 
at a distance of 140 p c. In comparis on, we show in Fig. [15] the 
v=l-0 CO line fluxes dBrittain et alj |2003) for the massive disc 
AB Aur and the low-mass disc HD141569A. The observed fluxes 
from HD141569A should b e scaled from its distance of 99 pc 
dvan den Ancker et alJI 19971) to 140 pc. Fluxes for the hot lines 
from the HD141569A disc are ~ 10-"* W mr~ after scaling the 
fluxes to 140 pc. When UV pumping is included, the hot line 
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10 100 I 10 100 I HI 100 

r[AU] r[AU] r[AU| 



Fig. 11. Model lax structures and abundances relative to H-nuclei for atomic He, electron, H, H2, H2exc (v >0), C^, C, CO, CH4, 
OH, and H2O. The contour plots in the ionized carbon abundance plot show log(^/n), where x is the enhancement with respect to 
the interstellar UV and n the total number of H-nuclei. In a few panels, dust temperature (Td), gas temperature (water abundance 
panel), or extinction Ay contours are overplotted by dashed lines. 



fluxes from our models are of the order of 10"'^-10"'^ W m"-. 
The shapes of the SLEDs for models with a small inner radius 
(7?in= 1 AU, models 1 & 3) reflect the effects of highly optically 
thick lines. On the other hand, the shapes of SLEDS of mod- 



els with a large inner radius (Rin- 20 AU, models 2 and 4) are 
typical of optically thin or moderately optically thick lines. The 
observed line fluxes are best matched by models with a small in- 
ner radius whereas the shapes of the observed SLED are closer to 
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Fig. 12. Location of the CO v = 1 - / = 1 9 - 1 8 P( 1 9) at 4. 844 |/m and CO v = 4 - 3 / = 20 - 1 9 P(20) at 4.692 yum line emissions 
for the Mdisc=10"^ Mq, /?in=l AU disc models with (left panels, model la) and without (right panels, model lb) UV-pumping. The 
solid black contours in the CO density panels encompass the regions that emit 49% of the fluxes. The black dashed-line contours 
contain 70% of the fluxes in the vertical direction. The white contours are at gas temperatures of 500 and lOOOK. Each panel is 
divided into three sub-panels. The upper sub-panels show in color scale the atomic hydrogen density structure, the middle-panels 
show the molecular hydrogen density structure, and the lower-panels show the electron density structure. 



models with a large inner radius. To obtain fluxes that are close 
to the observed values, the emitting gas should be close to the 
star, but at the same time lower surface densities are needed to 
have moderate optical depths for the CO lines. Interestingly, UV 
pumping will increase the line v'=l-0 fluxes only by a factor of 
few and the \'-4-3 fluxes by orders of magnitude. In our mod- 
els, UV pumping does not dramatically change the shape of the 
SLEDS for v =1-0 and v=4-3 transitions. 

We ran models with inclination from 1 to 89 degrees with 
the line ray-tracing module. The line fluxes vary by up to 30% 
for inclinations below 70 degrees. At high inclination, the cold 



CO in the outer flaring disc starts to reabsorb the emission from 
the inner disc. 



3.3.3. CO ro-vibrational level excitation 

The CO vibrational levels in the ground electronic state can be 
populated by collisions with the main chemical species (H, H2, 
electrons. He), by formation pumping (i.e. formation of CO in an 
vibrationally excited level), by pumping by dust emission, or by 
UV-fluorescence pumping. Th e region of high warm CO abun- 
dance is H rich and H2 poor dKamp & Dullemondll2004l) (see 
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Fig. 13. Similar to Fig.Obut for the Mdisc = 10"'* Mq, /?in=l AU disc models (3a and 3b). 



Fig.[T2land[T3]). We also plotted the excitation temperatures. We 
chose two typical transitions, one from the v = 1 level and one 
from the v - 4 level, to illustrate the efficiency of the excitation 
mechanisms. The relative population between the v = 1,7= 19 
and V = 1, / = 18 and between v = 4, 7 = 20 and v - 4,J - 19 
levels are shown as excitation temperatures in Fig. [17] and [18] for 
the high- and low-mass disc models, respectively. The excitation 
temperatures between the v = 1,/ = 19 and v = 0, 7 = 18 
levels in the line-emitting region are below the gas kinetic tem- 
peratures but higher than the radiation temperature at the line 
frequency 4.844 micron (subthermal population, see Fig. IA.2I 
IA.4l and lA.5l in the appendix). However, the excitation temper- 
atures between the v = 4, 7 = 20 and v = 3, 7 = 19 levels are 
much higher than the gas kinetic and radiation (at 4.692 micron) 
temperatures (suprathermal population). We can distinguish be- 



tween the excitation mechanisms of the v = 1 level and that of 
the I' > 1 levels. 

The V = 1 level is mostly populated by collisions and IR 
pumping because the disc densities in those emitting regions are 
close to the critical densities (lO^'-lO'^cm^^*). The v = 1 level 
population is weakly affected by UV-fluorescence pumping (see 
the upper panels of Fig.[T7land[T8ll. Likewise the line fluxes for 
the fundamental transitions change by at most a factor four when 
UV-fluorescence pumping is taken into account (Fig.[T4ll. 

The collision rates with electrons are ~ 100 times higher than 
with H but the electron abundance is ~ 10"^ lower than the H- 
atom abundance. The collision rates with He are of the order 
of 10"'^ cm-' s"' at 500 K. The abundance is 0.075 times less 
than H+H2. Therefore He is not a major collision partner. The 
main collision partner in the CO ro-vibrational emitting region 
is atomic hydrogen, whose de-excitation collision rates with CO 
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Fig. 14. CO V = 1 - continuum-subtracted line fluxes. The blue triangles show the fluxes for models with UV-pumping, whereas 
the red diamonds show the fluxes for models without UV-pumping. 
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Fig. 15. CO SLED from AB Aur (left panel) and HP 141569A discs (with the fluxes scaled to a distance of 140 pc, right panel). 
Both observations are taken from lBrittain et alj (|2003|) . 



are two orders of magnitude (10 cm^ s ' at 500 K) higher 
than rates with H2. 

The v > 1 levels have higher critical densities than the v = 1 
levels and require very high gas temperatures Tgas >1000 K and 
high densities to be populated efficiently by collisions. Pumping 



of the V > 1 levels by UV fluorescence becomes the promi- 
nent excitation mechanism. The excitation temperatures of UV- 
pumped levels (lower left-hand panels in Fig. [TTl andfTSl) tend 
towards the effective radiation temperature (brightness temper- 
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Fig. 16. CO V = 4 - 3 continuum-subtracted line fluxes. The blue triangles show the fluxes for models with UV-pumping whereas 
the red diamonds show the fluxes for models without UV-pumping. 



ature) in the UV (Fig. IA.31 the radiation temperatures at 0.15 
micron). 

The effect of UV pumping becomes more pronounced in the 
10""^ M0 disc models (3a and 3b, FigfT3]l. where the gas densi- 
ties and hence the UV-pumped high-v level quenching efficien- 
cies are lower and the UV photons are less absorbed by the dust 
grains (Fig. IA.3I ). In the Mdisc=10"^ Mq disc model, the UV 
pumping is confined to a very thin layer at the inner disc rim, 
whereas in the low-mass disc models (Mdisc=10"'* Mq), the UV 
pumping occurs till ~ 0.5 AU inside the disc because the UV 
photons can penetrate deeper in the disc. 

The UV-fluorescence pumping can affect dramatically the 
line fluxes (Fig.fTSIl. In particular the hot line v = 4 - 3 P(20) is 
two orders of magnitude stronger with UV-pumping. The effect 
of UV pumping is better captured by computing the vibrational 
temperature (see Sect. l3.3.5T l. The relative rotational level popu- 
lations within a vibrational level are less affected by UV pump- 
ing than the relative vibrational level populations (see the shape 
of the SLED in Fig. [Hi. 

An alternative excitation mechanism is the formation pump- 
ing of CO. The CO formation reaction via C + OH has an 
exothermicity of -6. 5 eV and a rate of fcchem ~10"'" cm-' s"' 
dZanchet et al.ll2007h . The excess energy is used to pump CO 
to high vibrational levels with a maximum population at v=10 



dBulut et al.ll201 ih . The formation of CO constitutes a pumping 
mechanism for CO. 

We can compare the efficiency of the chemical pumping with 
the efficiency of the population of the v=l level by collisions 
with atomic hydrogen by noting that OH reaches a maximum 
abundance of ~10"^ (Fig. fTTI). The chemical pumping rate as- 
suming atomic carbon abundance [C]~10""* is A:chem«c«OH ~ 
10""" «~ cm-' s' , where n is the gas density. Collision rates with 
atomic hydrogen are of the order of 10^'^n^ cm^ s', a few orders 
of magnitude more efficient than the chemical pumping. 

3.3.4. Rotational diagram for tine fundamental transitions 

The analysis of the line fluxes using a rotational diagram can pro- 
vide more insight than the study of the shape of the SLEDs. We 
plotted in Fig.[T9]the fundamental transition rotational diagram, 
i.e. the log of the level column density versus BcoJ'iJ' + 1), 
where Z?co is the rotational constant in Kelvins, for the 10"^ 
and 10 "* Mq /?in=l AU and 20 AU models with (blue symbols) 
and without UV pumping (red symbols). The shape of a rota- 
tional diagram is affected by NLTE (subthermal population and 
IR/UV pumping) and optical depth effects. To obtain the opti- 
cally thin emission in LTE (purple solid lines), we populated the 
CO ro-vibrational levels assuming Tex = Tgas- The line fluxes 
were computed in the optically thin approximation, i.e. optically 
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Fig. 17. . Excitation temperatures between the upper and lower level of the CO v = 1-0 P(19) transition at 4.844 jim v = 1,7 = 
19 (upper two panels) and that of the CO v = 4 - 3 P(2Q) transition at 4.692 fim v - 4-, J - 20 (lower two panels) transitions for the 
^disc=10"~ Mq, /?in=l AU disc models with (left panels, model la) and without (right panels, model lb) UV-pumping. The solid 
contours (in white or black) in the CO density panels encompass the regions that emit 49% of the fluxes. The black dashed-line 
contours contain 70% of the fluxes in the vertical direction. The white contours are at gas temperatures of 500 and lOOOK. 



depth effects, which would have prevented the flux from increas- 
ing indefinitely, were not taken into account. 

The optically thin CO LTE populations probe the actual ki- 
netic temperatures in the discs. The populations were calculated 
by summing for each level all the CO molecules at that level 
in the models, assuming that the population is in LTE. For each 
model, the rotational population of the v = 1 level can be rela- 
tively easily matched by a single temperature T^'^™ because the 
CO fundamental lines probe a small region in discs where the 
gas is warm and dense. Moreover, the densities are high enough 
for the low- and medium-7 rotational levels within the v = 1 
level to be close in rotational LTE. The line profiles will be wider 



for high-y than for low-7 transitions in discs not seen edge-on. 
The LTE temperatures (fits to the purple lines shown with black 
crosses) Tlte are given in Table H] The gas is warmer for lower 
mass, hence lower density discs because the molecules, which 
are efficient gas coolants, are more rapidly formed in dense 
gases. 

Since the stellar radiation decreases with distance squared, 
it is normal that the gas is cooler for discs with an inner radius 
starting at 20 AU rather than at 1 AU. Since CO ro-vibrational 
lines are efficient coolant for the gas, the gas probed by the CO 
gas is warmer (i.e. with higher T^^J^^) when the number of ro- 
vibrational levels decreases (Table SJi. 
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Fig. 18. . Similar to Fig.Obut for the Mdisc^lO""* Mq, R-.^^l AU disc models (3a and 3b). 



The red and blue lines in Fig. [191 correspond to the CO ro- 
vibrational NLTE level population computed locally in the disc 
using the 1 + lD escape probability method. Those lines show 
that the level populations are rotationally subthermal. The effects 
are stronger for less massive (10^^ Mq), hence less dense discs. 

The ro- vibrational NLTE population is orders of magnitude 
smaller than the ro-vibrational LTE population for all models, 
pointing to a subthermal population of the excited ro-vibrational 
levels; the rotational level populations within a vibrational level 
can be in LTE, while the vibrational level populations are sub- 
thermal. Assuming a transition probability of 10 s"' and a rate 
coefficient of lO"'*" cm^ s"' (H2) to 10"'^ cm^ s"' (H), the crit- 
ical density is 10"-10'^ cm"^. Such high densities are reached 
only close to the star in the inner disc midplane of massive discs. 

Radiative transfer of the optically thick lines decreases fur- 
ther the populations derived from the analysis of rotational dia- 



grams. One consequence is that rotational temperatures derived 
from fitting points from low-/ optically thick lines in rotational 
diagrams would appear lower than the actual gas kinetic temper- 
atures. In appendix lA.il we provide a simple analytical method 
to analyse observations of optically thick '^CO v=l-0 emission 
lines. Optical thickness results in flux differences between lines 
of transitions that have the same upper level but lower levels 
that differ by A/ and different transition probabilities: the P 
(J" - J' + I) and the R (/" - J' - I) branches. As expected, 
the column density differences between LTE and NLTE popula- 
tion are strong in models la and lb and are weak in model 4a 
and 4b. 

The /^-branch transitions generally have higher transition 
probabilities (i.e. with higher Einstein-A coefficients) than P- 
branch transitions. Hence, if the lines are optically thick, the 
column densities derived from /^-branch lines are the smallest 
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Fig. 19. Disc model CO v = 1 population diagrams: LTE and NLTE with and without UV pumping in the optically thin case (solid 
lines) and optically thick case (symbols). All the data have been normalized to the 7=0 population in LTE. The blue diamonds 
correspond to the column densities after line transfer in models with UV pumping whereas the red triangles correspond to column 
densities in line-transfer model results without UV pumping. The LTE populations are fitted by straight lines (black plus signs). The 
separations between P- and R- branches are clearly seen in the NLTE optically thick cases. 



(Nj' oc Fy/Aui and Fy are almost the same). Lines are optically 
thick up to very high J, giving rise to a steep slope for the low 
J and to a shallow slope for the moderate J lines (symbols). 
Observ ed rotational diagrams of Herbig Ae discs show similar 
shape s dBlake & Boogerj2004trBrittain et al.ll2007l:lBrown et alj 
I2OI2I) . Only hnes from J' > 40 (> 25 for the 10 Mq, 
^in=20 AU disc model) are optically thin: their populations join 
the optically thin NLTE populations and the differences between 
the P and R branch column densities fade away. The optical 
thickness of the lines and rotationally subthermal population of 
the levels can cause the lines to be actually optically thick but ef- 
fectively thin, especially in discs where the continuum emission 
at 4-5 micron is optically thin. 

UV pumping affects the level population for all levels with 
V > 1 whatever the inner disc radius. The main effect of UV 
pumping is to populate the high-7 and v > 1 (see discussion 
below and Fig. l20t levels at the expense of the lower v = 0, 7 
levels. We also ran the models with UV pumping but without 
collision excitations between the ground and electronic excited 
states. However, we noticed no difference over a few percents. 

The shape of the rotational diagrams (symbols) shows that 
it would be difficult to fit straight lines to estimate the excita- 



tion temperatures. Moreover, they do not reflect the actual CO 
population due to strong optical depth effects. 

3.3.5. Vibrational temperature and UV/IR fluorescence 
pumping 

The vibrational temperature measures the relative population of 
the vibrational levels. We plotted in Fig. [20l and 1211 the log of the 
vibrational level population (ln(A^,,)) as function of the vibration 
energy E^. = 3122v, where v is the vibrational quantum number 
(iBrittain et al.ll2007l) . The vibrational level population is the sum 
of all the rotational level populations at that vibrational level. 

We derived the vibrational temperatures by fitting straight 
lines through the data in models with and without UV pump- 
ing (nine and five vibrational levels) apart from the 10"^ Mq 
Rin-l AU model where only v=4, 5, and 6 levels are taken into 
account because a straight line could not be fitted when those 
levels were included. The values are reported in the last column 
of Table m 

Even without UV pumping, the vibrational temperatures are 
higher than the LTE rotational temperatures. However, the rota- 
tional temperatures derived using points in the shallow slope of 
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Fig. 20. CO vibrational diagram with and without UV pumping for the four models. We used the v-l-6 population after line 
radiative transfer. The blue diamonds correspond to models with UV pumping whereas the red triangles correspond to models 
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7?d,sc= 20 AU). 



the diagrams will yield temperatures higher than the UTE values 
and closer to the non-pumped vibrational temperatures. 

Two clear trends can be seen from the vibrational diagrams 
and the temperatures in Table |4] First, the vibrational tempera- 
ture is higher for the lower mass disc. Second, the vibrational 
temperature increases with the inner hole size if UV pumping 



is switched on, whereas it decreases if UV pumping is switched 
oflF. 

When the UV pumping is not present, the high vibrational 
levels can be populated by collisions and/or by IR pumping. The 
collision rates for vibrational transitions are lower than for rota- 
tional transitions. If the vibrational levels are populated by col- 
lision only, then the vibrational temperatures should be lower 
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Fig. 21. Similar to Fig.|20]but for models 5a, 5b to 8a, 8b. Among the models shown in this figure, the vibrational temperature 
differences between models with and without UV-pumping are the largest in the lowest mass disc with a inner hole (upper-right 
panel, model 6a with Mdisc - 10"^ Mq, Rdisc= 5 AU). 



then the rotational temperatures. Therefore, collisions play a mi- 
nor role in populating the high-v ro-vibrational levels, leaving IR 
pumping as the dominating mechanism. In the absence of UV 
pumping, the vibrational temperatures are still higher than the 
maximum dust temperatures in the discs (max(rdust) in Table|4]i. 
The IR continuum flux between 1 and 5 /im is dominated by 
the dust emission for the models with Rin-l AU (models 1 and 
3) and by the stellar emission for models with an inner hole 
(7?in=20 AU; models 2 and 4). The vibrational temperatures are 
the same whether five or nine vibrational levels are taken into 
account in the models, suggesting that the eff'ect of IR pumping 
does not depend on the vibrational level. 

As the gas cools with increasing inner radius, the high-v lev- 
els are less populated. Absorption of UV photons directly by 



CO molecules will favour the population of high-v levels unless 
high gas densities quench the overpopulation. The gas density 
at the inner edge decreases with the edge's distance to the star. 
Higher Rin values and lower mass discs will have less collisional 
de-excitations of UV-pumped high-i' levels. 

3.3.6. Effect of IR pumping on the rotational population of the 
V =1 level 

The left-hand panels in Fig.|22]show the differences between the 
case with two and five vibrational levels. For the massive disc 
with Mdisc =10"^ Mq, the population decreases as more vibra- 
tional levels are included. The CO v >1 levels can be pumped 
by absorptions in the hot and overtone bands of stellar or dust 
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Fig. 22. Disc CO v = rotational population at LTE and at NLTE for the R\n-i AU models with only the v = 0, v = 1,2, and 
V - 1 , 2, 3, 4, 5 levels taken into account. All populations have been normalized. 



IR photons, resulting in lower populations of the v=l, J lev- 
els (upper right-hand panel). The pure rotational levels in LTE 
cannot be explained by a single excitation temperature contrary 
to the rotational levels in the first vibrational excited state. The 
pure rotational transitions probe much larger disc areas with 
wide-ranging physical conditions than do the fundamental ro- 
vibrational transitions. This stems from the low critical densities 
needed to populate the pure-rotational levels. 

The right-hand panels in Fig. |22] compare the population in 
models where only rotational levels with v =0 are included and 
where the levels v =0 and v =1 are included. In the optically 
thin cases (solid lines), the differences between the cyan and blue 
lines reflect the rotationally subthermal populations for the high- 
J populations. The low-/ levels are at LTE. However, the high 
optical depths for the low- and mid-/ lines result in lower de- 
rived population and the appearance of a change in the slope at 
/=8-12. At high / {Erot >2000 K), the overpopulations when 
only the ground vibrational (v=0, /) are modelling artifacts be- 
cause overlapping energies occur between the v=0, high-/ level 
energies and that of the v=l, low-/ levels. For a gas at temper- 
ature r, both the high / and low / of the v =1 levels can and 
should be populated. In warm and hot gases, models that do in- 
clude the pure rotational levels only may overestimate their pop- 
ulations. 



4. Conclusion 

We have implemented a complete CO ro-vibrational molecular 
model in the ground and A'll electronic state in the radia- 
tive chemo-physical code ProDiMo. We gathered existing colli- 
sional rate coefficients and used scaling rules to extrapolate the 
missing ones. Collision rates of CO with atomic hydrogen are 
two orders of magnitude higher than rates with molecular hy- 
drogen and three orders of magnitude higher than with He. 

We ran models that include continuum, NLTE-line radiative 
transfer assuming a Voigt line profile, chemistry, and gas ther- 
mal balance to understand the effect of UV on the population 
of high-vibrational and high-rotational levels in protoplanetary 
discs. CO ro-vibrational lines are emitted within the first few 
AU and probe disc regions where CO is quickly formed via C + 
OH, with OH being formed via the fast reaction of vibrationally 
excited H2 with atomic oxygen. 

The amount of atomic hydrogen is high in the CO line- 
emitting region. The v = 1 level population is dominated by 
collisional excitation with atomic hydrogen with contributions 
from IR-pumping and UV-fluorescence pumping. Excitation of 
high-/ and high-\' levels by UV and IR photons can be efficient, 
especially in low mass discs where collisional quenching is less 
rapid than in high-mass discs. The molecular gas is rotationally 
"cool" and subthermal but vibrationally "hot" and suprathermal. 
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However, only the detection of hot lines (v' > 2, Av=l) in discs 
can provide a test of the efficiency of UV-fluorescence. 

The pure rotational transitions probe the entire disc while 
the ro-vibrational transitions probe the inner warm disc region. 
Therefore the rotational temperature depends on the vibrational 
number of the initial level. It is difficult to derive accurate tem- 
peratures and column densities from the optically thick '-CO 
fundamental emissions. Spatially and spectrally resolved CO, 
'^CO, C'^O, and C'^O ro-vibrational observations are needed to 
disentangle the optical depth effects and constrain the location 
of the line emissions. The low- and medium-/ rotational levels 
on the ground and first excited vibrational levels are populated 
in LTE in discs due to the low critical densities of the pure rota- 
tional transitions. 

Future studies will focus on studying the effects of disc pa- 
rameters (mass, size, flaring index, ...) on the CO ro-vibrational 
emissions and will include larger numbers of vibration levels and 
electronic excited states. 
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Appendix A: Appendix 

A.1. Rotational diagram analysis for levels in LIE and 
optically thick lines 

We adopted the rotational diagram analysis that takes optical 
dept h effe cts into account (also ca lled population diagram analy- 
sis injo^dsmith & Langerlll999h . A similar approach was taken 
by iGoto et al.l ( l201ll) . The line intensity emitted by a slab at a 
single temperature reads 



log ly = log (tvBv) + log 



1 



(A.l) 



By is the Planck function at the excitation Tex^Trot-Tyit- The 
second term in the right-hand side is an optical depth correction 
factor to the optically thin rotational diagram. The term /? = (!- 
exp(-Tv))/Tv is akin to the escape probability. The maximum 
optical depth occurs at Eroi ~ kTroi, where Eioi is the rotational 
upper level energy of an emission line. In the Born-Oppenheimer 
approximation, E,, = Erot + Eyn,. The rotational energy for CO 
can be approximated by Eioi - BroiJuiJi, + 1), with the rotational 
constant - 2.76 K. The maximum optical depth occurs when 
Ju - y/Trot/Brot- The line optical depth is 

" (^](^)('K^)^''^'''''^" " l)-^(^™"^'*)^(CO)/cos(0, 

(A.2) 

where v is the line frequency, A the Einstein spontaneous emis- 
sion probability of the transition, c the speed of light, Av the 
turbulent width, N{CO) the CO column density, / the disc incli- 
nation, and 

g-E,„,lkT„„ g-E,,i,/kT,,i, 

x{Tro„ T,ih) = — X — (A.3) 

'Jrot (Jvib 

is the CO fractional population in the initial level (v',7')- For 
CO, Qro, - kTrorlBco + 1/3 (Bco = 2.76 K) and Qy,,, = 
1/(1 - exp(-3122/r,.,fc) are the rotational and vibrational par- 
tition function respectively (Brittain et al. 2007). 

The rotational diagram of AB Aur (Fig. lA.ll ) shows three 
parts: a steep slope corresponding to increasing optical depth un- 
til Erot ~ kTex, then a shallow slope where the line optical depth 
decreases because the higher the rotational level, the less they 
are populated. At very high J, the slope steepens again. This be- 
haviour also appears in our theoretical rotational diagrams, but 
the second turning point occurs at higher J than in the observa- 
tions. The second slope change corresponds to lines with t <1. 
Assuming that the population is in LTE, we fitted the AB Aur 
rotational dia gram by a model that takes th e optical depth effects 
into account (iGoldsmith & Langeijll999h . The model parame- 
ters are the CO column density A^(CO)=4.2 x lO'^ cm"^, the 
mean excitation temperature Tex - 600 K (we assume that the 
vibrational temperature is equal to the rotational temperature), 
the turbulent width Av=0.05, and the inclination /=30 degree. 
In the upper panel of Fig. lA.ll the analytical solid curve com- 
pares well with the observations. The dashed-line curve shows 
the same model where optical depth effects are not taken into ac- 
count. The lower panel shows the derived optical depths, which 
reach ~ 47. 
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Fig. A.l. The upper panel shows the comparison between AB 
Aur and an analytical model rotational diagram for the '^CO 
y=l-0 transitions observed by Brittain et al. (2003). The solid 
blue line shows the fit by an analytical model that takes optical 
depths into account. The results from the same model but with- 
out the effect of optical depth are shown in dashed blue line. The 
lower panel shows the derived line optical depth. 

A.2. Excitation, radiation, and gas kinetic temperatures 

We define the radiation (brightness) temperature Trad at a given 
wavelength as the equivalent blackbody temperature that will 
match the specific intensity computed by the continuum radia- 
tive transfer at a given location in the disc (Fig. lA.2] and lA.3t . 
We also show temperature ratios in Fig. IA.4l and IA.5I the exci- 
tation over the radiation and the excitation over the gas kinetic 
temperatures. 
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Fig. A.2. Radiation temperature Trad at 0.15 micron and 4.844 micron for the Mdisc = 10"^ Mq, Rin-l AU disc models (left panel) 
and for the Mdisc - 10""* Mq, Rin-l AU disc models (right panel). The black contour shows the regions that emit 49% of the fluxes 
at 4.844 micron. The black dashed-line contours contain 70% of the fluxes in the vertical direction. 
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Fig. A.3. Radiation temperature Trad at 0.15 micron and 4.692 micron for the Mdisc = 10"^ M©, Rin-l AU disc models (left panel) 
and for the Mdisc - 10""* Mq, Rin-l AU disc models (right panel). The black contours show the regions that emit 49% of the fluxes 
at 4.692 micron. The black dashed-line contours contain 70% of the fluxes in the vertical direction. 
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Fig.A.4. TexCv = 1,/ = 19)/ri-ad(4.844) (upper panels) and TexCv = 1,7 = 19)/Tg^^ (lower panels) structures for the Mdisc=10"^ 
Mq, /?in=l AU disc models. The left panels correspond to the model with UV pumping and the right panels to models without UV 
pumping. The contours are labelled in red. The white contours encompass the regions that emit 49% of the fluxes at 4.844 micron. 
The dashed-line contours contain 70% of the fluxes in the vertical direction. 
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Fig.A.5. TexCv = 1,/ = 19)/ri-ad(4.844) (upper panels) and TexCv = 1,7 = 19)/Tg^^ (lower panels) structures for the Mdisc=10"'* 
Mq, /?in=l AU disc models. The left panels correspond to the model with UV pumping and the right panels to models without UV 
pumping. The contours are labelled in red. The white contours encompass the regions that emit 49% of the fluxes at 4.844 micron. 
The dashed-line contours contain 70% of the fluxes in the vertical direction. 



26 



